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Abstract 

To our best knowledge there is only one example of a lattice system with long-range 
two-body interactions whose ground states have been determined exactly: the one- 
dimensional lattice gas with purely repulsive and strictly convex interactions. Its 
ground-state particle configurations do not depend on the rate of decay of the inter- 
actions and are known as the generalized Wigner lattices or the most homogenenous 
particle configurations. The question of stability of this beautiful and universal re- 
sult against certain perturbations of the repulsive and convex interactions seems 
to be interesting by itself. Additional motivations for studying such perturbations 
come from surface physics (adsorbtion on crystal surfaces) and theories of corre- 
lated fermion systems (recent results on ground-state particle configurations of the 
one-dimensional spinless Falicov-Kimball model). As a first step we have studied a 
one-dimensional lattice gas whose two-body interactions are repulsive and strictly 
convex only from distance 2 on while its value at distance 1 is fixed near its value 
at infinity. We show that such a modification makes the ground-state particle con- 
figurations sensitive to the decay rate of the interactions: if it is fast enough, then 
particles form 2-particle lattice-connected aggregates that are distributed in the 
most homogeneous way. Consequently, despite breaking of the convexity property, 
the ground state exibits the feature known as the complete devil's staircase. 
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1 Introduction 



The present day equilibrium statistical mechanics is, without a great exaggeration, a 
theory of lattice systems with translation-invariant short-range interactions. The first step 
in a study of low-temperature properties of such systems amounts usually to determining 
their ground-state configurations [|l|]. Provided the underlying interactions are short- 
range, numerous methods of searching for ground-state configurations are available. In 
the case of one- dimensional lattice systems, there is even an algorithmic method [fj. |3|. 
In higher dimensions, let us mention the powerful method of m-potential which leads to 
success in many cases of interest |]. 

However, if the translation-invariant interactions are long-range ones, the situation 
is drastically different. The exact results are scarce. To our best knowledge they are 
available in two cases only: (1) a version of one-dimensional Frenkel-Kontorova model 
|4|] and (2) a one-dimensional lattice-gas model [[5], |J. This is in contrast with the case 
of continuous systems, where a number of results by various authors are available; an 
overview of these results and the corresponding references can be found in a review paper 
by Radin 0. 

We would like to mention that even these one-dimensional lattice cases are of great 
interest, being no academic problems but important scientific issues, well supported by 
the physics of quasi-one- dimensional materials highly anisotropic layered systems ||, 
and adsorbtion of molecules on crystal surfaces [II], [II], O. 

The lattice-gas models with purely repulsive and strictly convex two-body interactions 
emerged from considerations of orderings of electrons in quasi-one-dimensional conductors 
by Hubbard H, [5| and orderings of monolayers of atoms adsorbed on solids by Pokrovsky 
and Uimin ||7 The model still appears to be a cameo in this field. The periodic ground- 
state particle configurations of this model have been characterized exactly in || and ||. 
Particles are distributed as far as possible from each other, respecting restrictions imposed 
on their locations by the underlying lattice. These configurations, called by Hubbard the 
generalized Wigner lattices f|, are independent of any further details of the interaction 
potential. Moreover, the ground-state configurations exibit an interesting feature known 
as the complete devil's staircase |13|, [4| 

In view of such an impressing universality of the result obtained independently by 
Hubbard and by Pokrovsky and Uimin, the question of the stability of this result against 
some perturbations of the interactions seems natural and interesting from the point of 
view of statistical mechanics. Moreover, the interest in this question is supported by other 
domains of physics. In surface physics, a remarkable activity consists in studying, both 
experimentally and theoretically, orderings of molecules adsorbed on crystal surfaces. The 
phenomenon can be modelled by means of one- dimensional lattice gases [|II], [11], |12| with 
specific two-body interactions, which sometimes constitute certain local perturbations 
of strictly convex repulsive interactions [III. Another motivation for studying effects of 



modifications of repulsive and strictly convex two-body interactions stems from the re- 
cent results concerning ground states of a version of the Falicov-Kimball model [TJ| - the 
one-dimensional spinless Falicov-Kimball model. The model can be thought of as a model 
of quantum itinerant electrons and classical localized particles called f-electrons, nuclei 
15], [16| , with only on-site interactions whose strength is the unique parameter of 



or ions 



the model - a sort of the Ising model in the field of correlated fermion systems. Such a 
system can be transformed into a classical lattice-gas model with fairly complicated long- 
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range and many-body interactions |16|]. Many years after the discovery of the generalized 
Wigner lattices, Lemberger |T7j found the periodic ground-state configurations of the lo- 
calized particles in the large-coupling one-dimensional spinless Falicov-Kimball model. He 
named them the most homogeneous configurations because of the special role they play 
in his ingenious procedure of differentiation and integration of particle configurations. It 
turns out that the configurations found by Lemberger are just the generalized Wigner 
lattices. Despite the complicated nature of the effective interaction between electrons and 



ions, for some details see [18], there is a numerical evidence [fL9|, |20| (based on the approx- 
imate method of restricted phase diagrams) that in the strong-coupling regime it can be 
approximated by a repulsive and strictly convex two-body potential that has the same 
set of periodic ground-state configurations. Further studies of the ground-state phase 
diagram of the one-dimensional spinless Falicov-Kimball model, carried out in [|l^, |20[| , 
revealed for medium and small couplings a number of new families of good candidates for 
ground-state configurations. Among them are the so-called n-molecule most homogeneous 



configurations, found for the first time and studied in detail in Roughly speaking 



they are obtained from the most homogeneous configurations by replacing single parti- 
cles by lattice-connected aggregates of n = 2, 3, . . . particles. The arguments in favor 
of the existence of n-molecule most homogeneous configurations, provided in [ I9| , are 
based not only on the approximate method of restricted phase diagrams but also on the 
study of the interaction energy between a few ions only. It has been found that for large 
values of the unique parameter, for which the most homogeneous configurations are the 
ground-state configurations, the interaction energy between a few, say two or three, ions 
is strictly convex. On the other hand, for those values of the unique parameter for which 
the n-molecule most homogeneous configurations are the ground-state configurations, the 
interaction energy between two or three ions gets considerably lowered at short distances 
(which apparently encourages forming n-molecules) and consequently it becomes noncon- 
vex at short distances. These results led us naturally to the question whether it is possible 
to obtain the n-molecule most homogeneous configurations as ground-state configurations 
of a lattice gas with two-body interactions, by modyfing, at short distances only, the two- 
body interactions whose ground-state configurations are the most homogeneous ones. 

The answer to this question is in the affirmative, at least in the case n = 2 which 
we have considered as a first step of our investigations. Namely, we prove that repulsive 
interactions whose value at distance 1 is set near their value at infinity, that are strictly 
convex from distance 2 on, and whose rate of decay is fast enough (according to a simple 
criterion), do the job. 

Needless to say that we believe the sort of result we obtained to be valid for any n, 
however a proof along the lines of the case n = 2 seems to be technically complicated. 

The paper is organized as follows. In Section 2, we introduce the system under consid- 
eration and give basic definitions. Then, in Section 3 we formulate our hypothesis, provide 
examples that reveal problems that have to be solved on the way towards the final result, 
and draw our strategy of proving the hypothesis. A relatively simple part of this strategy, 
Lemma 0, is proved in this section. After that we prove our main theorem, Theorem 1, 
by means of Lemma 0, 1, and 2. The proofs of Lemma 1 and 2 are given in Section 4 and 
constitute the major part of the paper. Finally, in Section 5 we provide a discussion of 
the obtained results, limitations and possible extensions. Some technical definitions and 
statements that are used throughout the text are collected in the Appendix. 
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2 Lattice-gas models and devil's staircases 



A classical lattice-gas model, considered here, is a system in which every site of a one- 
dimensional lattice Z can be occupied by one particle or be empty. Then, an infinite-lattice 
configuration is an assignment of particles to lattice sites, i.e., an element of Q — {1, 0} z . 
If X G £1 and A C Z, then we denote by X\ a restriction of X to A. We assume that 
the particles interact only through two-body forces and to a pair of particles at lattice 
sites i and j, whose distance is \i — j\, we assign the translation-invariant interaction 
energy V(\i — The corresponding two-body potential reads: V(\i — j)\)si(X)sj(X), 
where {si(X),i G Z} are occupations of sites in a configuration X; S{(X) assumes value 
1 if in the configuration X the site % is occupied by a particle and otherwise value 0. In 
terms of the above defined potential, the Hamiltonian of our system in a bounded region 
A amounts to the sum of the potential over all pairs of sites having nonvoid intersection 
with A: 

H A (X)= Y, n\i-3\)si(X) Sj (X). (1) 

{{i,jHM'}nA^0} 

For Y,X G Q, we say that Y is a local excitation of X, and write Y ~ X, if there 
exists a bounded AcZ such that X = Y outside A. 
For Y ~ X, the relative Hamiltonian is defined by 

H(Y,X) = Yl (y(\i-j\>i(Y) Sj (Y) - V(\i-j\) Si (X) Sj (X)) . (2) 

{{i,i}:{iJ}(W0} 

We say that X G £1 is a ground-state configuration of H if 

H{Y,X) > for any Y ~ X. 

That is, we cannot lower the energy of a ground-state configuration by changing it 
locally. 

The energy density e(X) of a configuration X is 

where |A| is the number of lattice sites in A. In a similar way we define the particle density 
in a configuration X, denoted p{X), simply the Hamiltonian in (^) has to be replaced by 
the number of particles in the region A. 

It can be shown that any ground-state configuration has the minimal energy density 
which means that local conditions present in the definition of a ground-state configuration 
enforce global minimum of the energy density |2~T| 



For certain values of external parameters, like the chemical potential or the parti- 
cle density (depending on the used Gibbs ensemble), our models do not have periodic 
ground-state configurations. However, for any fixed value of such an external parame- 
ter, all ground-state configurations belong to one local isomorphism class. It means that 
they cannot be locally distinguished one from another. Every local pattern of particles 
present in one ground-state configuration appears in any other within a bounded distance. 
More formally, there exists a unique translation-invariant probability measure supported 
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by ground-state configurations. It is then necessarily the zero-temperature limit of equi- 
librium states (i.e., translation-invariant Gibbs states) |22f . 

In this paper, by the ground state of a model we mean precisely the above defined 
probability measure. 

If a system has a unique periodic ground-state configuration and its translations (this 
happens in our models for those values of the chemical potential that fix a rational particle 
density), then a unique ground-state measure assigns an equal probability to all these 
translations. For example, the Ising antiferromagnet has two alternating ground-state 
configurations but only one ground-state measure which assigns probability 1/2 to both 
of them. 

In the nonperiodic probability ground-state measure P gives equal weights 

to all ground-state configurations and can be obtained as a limit of averaging over a 
given ground-state configuration X and its translations r a X by lattice vectors a G Z: 
P = limA^z m J2aeA $( T a.X), where 5{r a X) is a probability measure assigning probability 

1 tO T a X. 

One more remark concerning ground states discussed here is in order. We consider 



exclusively the ground-state measures that are strictly ergodic [£3[ |J . In particular, every 



ground-state configuration in their support has uniformly defined densities of all local 
particle configurations. Moreover, if a local particle configuration occurs, then it occurs 
with a positive density. That is to say, we do not consider ground-state configurations 
with interfaces (like kink ground-state configurations in the Ising model). 

We say that a set A C Z is lattice- connected if for every pair of lattice sites i, j G A 
there is a sequence i\, . . . , i n such that i\ — i, i n — j and ik, ik+i, k — 1, . . . , n— 1 are the 
nearest-neighbor sites. 

In the sequel, in order to describe some local configurations X\, Ya, for a lattice- 
connected A, we find it convenient to introduce on Z a coordinate-axis with lattice sites 
located at integer points and to call the positive direction the right one while the opposite 
direction - the left one. The environment of A, Z — A, splits naturally into the left 
environment, consisting of the sites preceding A, and the right environment, consisting 
of the sites following A. It will also be convenient to set the zero of the axis at the last 
occupied site of the left environment so that the position of a particle in A is positive and 
coincides with its distance to the left environment. The coordinate-axis introduced will 
be briefly called x-axis. 

We assume that the interaction energy of two particles separated by distance r, V(r), 
is summable and for r > ro it is positive and strictly convex (see the Appendix), hence 
decreasing. Such an interaction energy is denoted by V ro (r). It is convenient to normalize 
V ro {r) in such a way that lim^oo V ro (r) = 0. 

Then, we define the interaction energy of a particle at a site r in A with the left 
environment of A, V^(r), as the sum of interaction energies V ro (r — j) over all occupied 
sites j in the left environment. The properties of the interaction energy V rQ imply that, 
for r > ro, the function V r ^(r) is positive, strictly convex, and decreasing. 

In what follows some particle configurations play a distinguished part. Among local 
configurations that appear frequently in our considerations are atoms, i.e., occupied sites 
whose nearest-neighbor sites are empty. The location of an atom is identified with the 
location of the occupied site. Another local configuration is a lattice-connected set of n, 
n > 2, occupied sites whose left and right nearest-neighbor sites are empty. It is called 
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the n- molecule. The location of a n-molecule is identified with the location of the first 
particle of the molecule, i.e., the one with the smallest x-coordinate. Consequently, the 
distance between two n-molecules is identified with the distance between the first particle 
of one n-molecule and the first particle of the other n-molecule. 

However, our attention is focused on global particle configurations which can be char- 
acterized as follows. For every particle density p, there is a unique sequence of natural 
numbers d n , such that the separations between any pair of n-th nearest-neighbor particles 
are d n or d n + 1. Such configurations have been called the generalized Wigner lattices || 



or the most homogeneous configurations [17]. If the particle density is rational, then the 



corresponding most homogeneous configurations are periodic (the particle locations can 
be given by means of a construction given for instance in |J) while for irrational densities 
they are nonperiodic. 

At least in this paper, the main reason of interest in the most homogeneous configu- 
rations stems from the following theorem: 

Theorem (Most homogeneous ground-state configurations) 

In the canonical ensemble, i.e., for a given particle density p, the ground-state config- 
urations of a lattice gas (]]]) with an interaction energy V\ are the most homogeneous 
configurations. 

This statement has been proven (or at least a proof has been outlined) by Hubbard 
|| and Pokrovsky and Uimin ||. The ground-state phase diagram in the grand-canonical 
ensemble has been calculated heuristically by Bak amd Bruinsma [[R| while a proof for 



the related Frenkel-Kontorova model has been provided by Aubry []4[ and adapted to the 
lattice-gas model case by Miejdsz and Radin [^5|; it is outlined below. 

In the grand-canonical ensemble, to find the energy density of a ground state we have 
to minimize 

f{p) = e{p)-pp. (4) 

Now, e(p), i.e., the ground-state energy density for the particle density p, is differentiable 
at every irrational p and is nondifferentiable at any rational p 0, However, as a 
convex function, it has a left derivative d~e(p)/dp and a right derivative d + e(p)/dp at 
every p. It follows, that to have a ground state with an irrational density p of particles, 
one has to fix p(p) = de(p)/dp. For any rational p, there is a closed interval of chemical 
potentials [d~e(p)/dp,d + e(p)/dp], where the most homogeneous configurations of density 
p are the ground-state configurations. One can show that the sum of lengths of these 
intervals amounts to the length of the interval begining at the end of the half-line, where 
the vacuum is the only ground-state configuration, and ending at the begining of the 
half-line, where the completely filled configuration is the only ground-state configuration. 

As we have already mentioned, for any rational p, there is a unique (up to translations) 
periodic ground-state configuration with that density of particles - there is a unique 
ground-state measure. For any irrational p, there are uncountably many ground-state 



configurations which are the most homogeneous configurations. It has been shown in [[26 
that there is still the unique ground-state measure supported by them. 

The particle density versus the chemical potential of particles, p(p), is constant in each 
set of this partition. Moreover, it is a continuous function on the real line and is inversion 
symmetric with respect to the point (po, p{po)), where /iq is the chemical potential for 
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which the free energy density is hole-particle symmetric (p(//o) — 1/2)- The curve is 
classified as a fractal one and named the complete devil's staircase 

The last remark is that without any loss of generality we can restrict our considerations 
to systems whose particle density p does not exceed 1/2. Then the most homogeneous 
configurations consist exclusively of atoms. 



3 Basic ideas and the main result 

Whether we follow the argument of Hubbard [[J (based on a version of our Lemma Al) or 
the argument of Pokrovsky and Uimin || that proves Theorem 0, we find that such a proof 
consists essentially of two stages. In the first stage we "chop" configurations that contain 
n-molecules with some n = 2, 3, ... off the set of all configurations with particle density 
p, so we are left only with configurations that are composed of atoms whose density is p. 
The second stage is like a "fine tuning" that amounts to precise adjusting the positions 
of atoms in order to minimize the energy density. 

Is the strategy outlined above useful if the two-body strictly convex for r > 1 interac- 
tion energy V\, that appears in Theorem 0, is replaced by an interaction energy V ro which 
is strictly convex only for r > r , r = 2, 3, . . .? 

Consider first the second stage. Suppose that the set of configurations composed of 
atoms only is replaced by the set of configurations that consist exclusively of n-molecules 
with fixed n, separated by at least distance r Q + n — 1. Let the particle density of these 
n-molecule configurations be p. This class of configurations we denote by C r n o p . Clearly, 
the problem of determining the ground-state configurations in C™ can be reduced to the 
analogous problem but in the class C* o , n and with V ro replaced by the effective two-body 

interaction between n-molecules, Vff . This effective interaction can be naturally defined 
as the sum of interactions V ro between ordered pairs of particles such that the first member 
of a pair belongs to one n-molecule while the second member - to the other n-molecule: 

V r f(r) = nV ro (r) + (n - 1) (V ro (r + 1) + V ro (r - 1)) + 
(n - 2) (V ro (r + 2) + V ro (r - 2)) + . . . + {V ro (r - n + 1) + V ro (r +n - 1)) , (5) 

where r stands for the distance between the two considered n-molecules. Since the func- 
tion V^p(r) is a sum of functions that are strictly convex for r > r , it is also strictly 
convex for r > r . This implies that in the ground-state configurations, the first particles 
of n-molecules form the most homogeneous configurations with particle density p/n, and 
the same holds true for the second particles of n-molecules, etc. Thus, it is natural to call 
the obtained ground-state configurations the n-molecule most homogeneous configurations 
of particle density p. We summarize the above considerations in the following lemma: 

Lemma (Ground-state configurations restricted to C^ Q p ) 

In a lattice gas (|l|) with an interaction energy V ro , the ground-state configurations re- 
stricted to C,™ are the n-molecule most homogeneous configurations. 

Having generalized succesfully the stage two, it is tempting to turn to the stage one. 
Can we modify V\ in such a way that the ground-state configurations are in ? The 



general suggestion that comes from studies of the Falicov-Kimball model |19| is to set 
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the values of the two-body interaction energy at short distances close to zero. One might 
expect that setting Vi(l) = will force the system to form 2-molecules exclusively in the 
ground-state configurations. However, in order to arrive at such a result, one has to deal 
with other features of the interaction energy. At this point it is instructive to turn to 
examples. 

Example 1 

Consider the convex two-body interaction energy V that in some units is given by 1/(1) = 

0, V(2) = 4, V(3) = 2, 1/(4) = 1, V(r) = 0, for r > 5, and two periodic configurations 
(period 17) of particle density p = 8/17 whose elementary cells are of the form: [• • o oo • 
mooo»mooo]! and [•ioomoomoomooo] 2) where • stands for a particle while o - 
for an empty site. The energy per elementary cell in the first case is 21/(2) + 31/(4) = 11, 
while in the second case it is 31/(3) + 71/(4) = 13. Thus the configuration [.] 2 , which 
consists exclusively of 2-molecules, looses against the configuration which consists of 
2-molecules and 3-molecules. While the interaction energy chosen is not strictly convex, 
it is easy to see that the above result remains true if we do not set V(r) = for r > 5 
but extend V(r) to some Vi, with V\{r) positive and arbitrarily small from r = 6 on. We 
deal with the above problem in Lemma 2 which tells us that to exclude n-molecules with 
n > 3 from competition one should impose a condition on the relative strength of 14(2) 
with respect to V\(r) with r > 3. Namely, it is sufficient to require that the energy of the 
"tail" of V 2 , defined as W = Y^L^V 2 {r), is weak enough compared to V 2 (2) + V 2 (l)/2, 

1. e., V 2 {2) + V 2 (l)/2 > 7W/2. 

Example 2 

This example shows that configurations containing atoms are unlikely to be the ground- 
state ones. Let the two periodic configurations of particle density p = 3/11 be of the 
form: [• • (5o) • «(5o) • •(6o)] 3 (period 22) and [• • (4o) • (4o)] 4 (period 11). The notation 
"(5o)" stands for five empty sites in a row, etc. If V(r) = for r > 8, then the energy 
per cell of 22 sites amounts to 21/(6) + 51/(7) + 41/(8), in the case of [.] 3 , and to 41/(5) 
+ 41/(6), in the case of [.] 4 . Now, let the interaction energy V 2 be chosen to coincide with 
the finite-range interaction energy defined as follows: 1/(1) is nonpositive, 1/(2) = 28, 
1/(3) = 21, 1/(4) = 15, 1/(5) = 10, 1/(6) = 6, 1/(7) = 3, 1/(8) = 1, in some units, while for 
r > 9, V 2 (r) is an arbitrarily small positive extension of this function to a strictly convex 
function. For such V 2 the configuration [.] 3 wins overwhelmingly. 

Summing up, in order to implement the stage two in the case of a nonconvex interac- 
tion energy V 2 , we have to force the system, by modyfying its interactions, to make the 
configurations containing n-molecules with n > 3,4,... unfavourable energetically, and 
to show generally that for such interactions the configurations containing atoms cannot 
be among the ground-state ones. The latter problem is dealt with in Lemma 1. 

Lemma 1 applied for instance to the configuration [J4 tells us that, for any V 2 such 
that V 2 (l) < 0, if we take twice as large elemetary cell and rearrange the two atoms and 
the 2-molecule separating them, so that they form two 2-molecules distributed as follows: 
[• • (60) • «(4o) • •(6o)] 5 , then we win the energy. By Lemma 0, we can win even some 
more energy by adjusting the distances between the 2-molecules so that the resulting 
configuration is the 2-molecule most homogeneous one: [• • (5o) • «(6o) • •(5o)] 6 . 

What we have said above, summarized in Lemma 0,1, and 2, shows that the strategy 
that lead to the proof of Theorem 0, which was concerned with strictly convex interaction 
energies, can be applied also to some nonconvex interaction energies V 2 . 
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Theorem 1 (2-molecule most homogeneous ground-state configurations) 
Consider a lattice gas ([[]) with a nonconvex interaction energy Vi and a particle density 
p < 1/2. If V 2 (l) = and V 2 (2) > 7W/2, then the ground-state configurations are the 
2-molecule most homogeneous configurations of particle density p. 

Corollary (Complete devil's staircase) 

In a lattice gas (0) whose interaction energy Vi satisfies the conditions given in Theorem 
1, the particle density versus the chemical potential of particles, p(/i), exibits the complete 
devil's-staircase structure. 

4 Two lemmas on eliminating atoms and n-molecules 
with n > 3 

Lemma 1 (Eliminating of atoms) 

Consider a lattice gas (|l|) with an interaction energy V2, such that V^l) < 0. Let X be 
a configuration that does not contain n-molecules with n > 3. Suppose that the local 
configuration Xa, where A is a bounded lattice-connected subset of Z, contains two atoms 
separated by k — 0, 1, 2, . . . molecules. Then, by rearranging the positions of particles in 
Xa, it is possible to construct a configuration Y such that p(Y) = p(X), Y"z-A = Xz_a, 
Fa consists exclusively of k + 1 2- molecules and e(Y) < e(X). 

Corollary: Among the configurations that do not contain n-molecules with n > 3, the 
lowest energy configuration consists exclusively of 2-molecules. 

Proof: We propose a proof by reductio ad absurdum. Let A be a bounded lattice- 
connected subset of Z and suppose that a local configuration Xa contains two atoms 
separated by k — 0, 1, 2, ... 2-molecules while Xz-a is arbitrary. The idea is to construct 
another configuration Y, with Yz-A = Xz-a and Ya consisting of k + 1 2-molecules 
exclusively, such that the total energy change, H(Y, X), is negative. 

We start with describing Xa- Going from left to right, its 1-st (or left) atom is 
separated by a empty sites from the left environment, after the first atom we meet the 

1- st 2-molecule separated from the 1-st atom by rii empty sites, later on - the second 

2- molecule separated from the 1-st one by n<i empty sites and so on. The (i — l)-th and 
the i-th 2-molecules are separated by n^ empty sites. Finally the 2-nd (or right) atom is 
separated from the last 2-molecule, i.e., the k-th 2-molecule, by n k+1 empty sites and by 
b empty sites from the right environment. 

As it has been said above, in the new configuration Y the configurations of the left 
and right environments remain unchanged while Ya is obtained as follows. We start with 
the 1-st atom and the left particle of the 1-st 2-molecule in Xa and put them at the two 
nearest-neighbor sites located in the center of the gap of rii empty sites. Clearly, if ri\ 
is even, then the position of the new 2-molecule (the 1-st 2-molecule in Ya)) is defined 
uniquely, the center of gravity of the 2-molecule coincides with the geometric center of 
the gap. On the other hand, if n\ is odd, then we can speak of the left and right central 
positions of the 2-molecule. In the first case the center of gravity of the 2-molecule is 
shifted by half the lattice constant to the left of the geometric center of the gap while in 
the second case - to the right. In the next step we create the 2-nd 2-molecule of Ya by 
taking the right particle of the 1-st 2-molecule and the left particle of the 2-nd 2-molecule 
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in X A and putting them in the center of the gap made by 712 empty sites. The procedure 
described is continued until k + 1 2-molecules in Y\ is created. Each time the gap between 
the 2-molecules in X A is odd we have a choice of two central positions. Therefore we may 
end up with as much as 2 k+1 different configurations. In the sequel we shall restrict 
ourselves to two simple choices, namely either we always choose the left central position 
and the corresponding Y is labeled Y L or we always choose the right central position what 
leads to the configuration Y R . 

Our aim is to estimate the energy changes related to passing from X to Y L , H(Y L , X) 
and from X to Yr, H(Yr, X). It is convenient to split the total energy of X A , E A (X), into 
two parts: E A (X) = E A nt (X) + E A :t (X). Here E A nt (X) is the internal energy, that is the 
sum of pair interaction energies V ro (r) over all pairs of particles in X\. The external energy 
of X A , E e A xt {X), can in turn be represented as Ef\X) = E%(X) + E%(X), where Efr{X) 
is the sum of pair interaction energies over all pairs of particles consisting of one particle 
in the left environment at configuration Xz~a and one particle in A at configuration X\, 
and E A (X) is defined similarly. Therefore 

H(Y,X) = E A (Y)-E A (X) = 
= E A (Y) - E L A {X) + E A (Y) - E*(X) + E A nt {Y) - E**(X), (6) 

for some Y. 

The remaining part of our proof consists of two stages. In the stage one we study only 
the external energy differences while in the stage two the internal energy differences. 

Stage one: estimating external energy variation 

Consider the pair of particles that constitute the i-th 2- molecule in Yl- Let E^(Yl) 
and Ej^(X) stand for the external interaction energy of that pair of particles with the left 
environment at configurations Yl and X, respectively. Then 

fc+i 

E L A (Y L ) - E L A (X) = £ (E^Y L ) - E?{X)) (7) 
i=i 

and similarly 

fc+i 

ERY L ) - E*(X) = E {E*{Y L ) - E*{X)) (8) 
i=i 

Therefore we are left with estimating the energy differences E^(Yl)—E^(X) and E^(Yl) — 
Ef'iyX) what requires a detailed description of positions, in configurations X and Yl, of 
the two particles that constitute the i-th 2-molecule in Y L . In X, the particles mentioned 
are the right particle of the (i — l)-th 2-molecule whose position according to the x- 
axis is denoted by Xj and the left particle of the i-th 2-molecule whose position is then 
Xi + ni + 1. Let us recall that the positions mentioned coincide with the distances between 
the particles considered and the left environment. In Yl the corresponding positions are 
Xi + [rii/2] and Xi + [rij/2] + I. Now we are ready to express the energies E[(Y L ) and 
E^(X) in terms of the interaction energy, V r ^(r), of a particle at r e A with the left 
environment (see Section 2 for the definition and properties). Then 

E^(X) = V L (x t ) + V L (x t + n t + l) (9) 
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and 

Ef{Y L ) = V L (xi + [m/2]) + V L { Xl + [m/2] + 1). (10) 

Therefore 

Et(Y L )-Et(X) = 

(v L ( Xl + [m/2]) - v L ( Xi )) + (y L ( Xi + [m/2] + 1) - v L ( Xl + m + 1)) . (li) 

Since V L (r) is a decreasing function of distance r, the inspection of the distances that 
appear in eq.(|TT|) shows that the first energy difference is negative while the second one is 
positive. Hence the net outcome can only be established by calling additional properties 
of V L {r). We use the convexity property of V L {r) (for r > 2). Estimating E^(Yl) — 
E^(X) amounts to estimating the variations of V L {r) at the intervals [xi,Xi+ [m/2]] 
and [ X i + [m/2] + l,Xi + m + 1] (see sec. 2). First suppose that m is even. Then both 
intervals are of the same length, m/2, but the second one is shifted to the right with 
respect to the first one. Therefore by convexity 

V L {xi + m/2 + 1) - V L { Xi + m + 1) < V L { Xi ) - V L { Xi + m/2) (12) 

and consequently 

E^(Y L ) - E^(X) < 0, (13) 

which is the desired result. However if m is odd, one of the intervals, the one located more 
to the right, is longer by 1 and an analogous reasoning does not reproduce inequality flT3]). 
Translating the longer interval to the left, so that its left end coincides with the left end 
of the shorter one, we obtain by convexity 

E?(Y L ) - E^(X) < V L ( Xi + [m/2]) - V L { Xl + {m/2}). (14) 

Remark: In a similar way one can derive a stronger inequality, namely the r.h.s. of ([14]) 
can be replaced by V( Xi + m) — V( Xi + m + 1). 

To get an estimate of the external energy E^iYi) — E e t ^ t {X) it remains to find the 
counterparts of inequalities ( [i~3|) and (|I4"D in the case of the right environment. For this 
purpose it is convenient to introduce a second coordinate-axis, the ?/-axis, similar to the 
x-axis but pointing in the opposite direction. We keep using the notions of left and 
right according to the a;-axis. The zero of the y-axis is set at the 1-st (i.e., most to the 
left) particle of the right environment, so that again the y-position of a particle in A 
coincides with the distance between that particle and the right environment. Now a while 
of reflection enables us to realize that the interaction energies with the right environment 
can be obtained from the interaction energies with the left environment by replacing Xi 
by yi and interchanging the square and curly brackets. Thus the counterparts of (Q) and 
(0) read 

E?{X) = V R { Vi ) + V R { Vi + m + l) (15) 

and 

E R {Y L ) = V R (y t + {m/2}) + V R (y t + {m/2} + 1), (16) 



11 



respectively, and the counterpart of ( |TT| ) is 

E?(Xl) - E*{X) = 

(V R ( Vi + K/2}) - V R { Vi )) + (V R ( yi + { ni /2} + 1) - V R (vi + n t + 1)) . (17) 
Again, if Hi is even we find 

E R (Y L ) - E R {X) < (18) 

and if rij is odd: 

E R (Y L ) - E R (X) < V R (y t + {n t /2}) - V R (y t + [ni/2]). (19) 

Summing up, the external energy variation on passing from X to Yl has an upper bound 
of the form 

Ef\Y L ) - ET\X) = E L A {Y L ) - E L A (X) + E R {Y L ) - E R (X) = 

fc+1 . fc+1 

< 

I 't \ J -' / I \ /I 

8=1 ' 1=1 

fc+1 



e - ^po) + e K(^) - *?po; 

8=1 8=1 

E (^(a* + [ni/2]) - K L (x, + K/2})) + (V R ( yi + {n,/2}) - V R { Vi + [ni/2])) . (20) 



8=1 



Again a while of reflection enables us to write down the upper bound for the change in 
the total external energy on passing form X to Yr, that represents a counterpart of (|20|) . 
It is enough to interchange in ( p0|) Xi with ?/j and R with L: 



e^(^) - TO = ^(y fl ) - ££P0 + <(f^) - <(x) = 
E Of(^) - ^PO) + E (^f(^) - £,*(*)) < 

8=1 8=1 



fc+1 



E (V R {vi + [ni/2]) - V*^ + {ni/2})) + (v L ( Xi + K/2}) - V\x i + [ni/2])) . (21) 



8 = 1 



By inspection of (E^) and (ETT) one finds that the upper bound given in (B^) is just the 
opposite of the upper bound given in (|2~lD . Thus, if .E^C^l) — E A xt (X) > 0, then its upper 
bound (p0|) is strictly greater than zero, consequently the upper bound fl2lD is strictly less 
than zero what implies that E^^Yr) — E A xt (X) < and vice versa. Therefore we proved 
that either E e A xt (Y L ) - E e A xt (X) or E e A xt (Y R ) - E e A xt (X) is nonpositive. 

Remark: Let us note here that the lemma is valid also if k = 0. 

Stage two: estimating internal energy variation 

Now we are going to estimate the internal energy difference E^^Yl) ~E l A t {X) ) which 
appears to be a considerably more laborious task than it was the case for the external 
energy difference. To arrive at a satisfactory upper bound of E A xt {YL) — E e A xt (X) we 
represented the external energy E A xt (Yi,) and E A xt (X) by a sum of external energies of 
pairs of particles in A that constitute 2- molecules in the configuration Yl (Yr). Such 
a partition of the external energy led to energy differences that could be estimated by 
means of the convexity property of the interaction involved. We shall need a counterpart 
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of this partition for E^*(Yj,) and E % ^ t (X). This time a natural object is not a pair of 
particles but two pairs of particles that constitute two 1-st nearest-neighbor 2-molecules 
in Yl (Yr), two 2-nd nearest-neighbor 2-molecules, and so on, until a pair of (k + l)-th 
nearest-neighbor 2-molecules. We can define in a natural way the internal interaction 
energy of those objects, in terms of V(r). First consider a pair of 1-st nearest neighbor 
2-molecules in Yl, say the i-th and (i + l)-th one, i — 2, . . . , k — 1 (so that the boundary 
molecules, the 1-st and the (k + l)-th one, are excluded from our consideration for a 
moment). To build up the i-th 2-molecule in Yl, we pick up the right particle of the 
(i — l)-th molecule in X and move it to the right by distance [rii/2], then the left particle 
of the i-th molecule in X and move it by distance {nj/2} to the left. Thus the i-th and 
(i + l)-th molecules in Y L are separated by the distance x^' + 1, where 

xW := g(rn,n».i) = {n t /2} + [n i+1 /2], (22) 

and we used the function g(s, t) defined in the Appendix. The internal interaction energy 
of the pair of molecules considered is defined as the sum of all the four interaction energies 
V(r) between the particles separated only by the gap a;- 1 ^: 

E?\Yl) = E YL ( ni ,n i+1 ) := (v(x^ + 1) + V{x { p + 2)) + (V(x^ + 2) + V{xf ] + 3)) . 

(23) 

While according to our description above the definitions fl2"2"|) and ( |2"3"D do not apply to 
these cases % = 1 and i = k, a direct inspection shows that 

= 0(711,712), 4 1} = g{n k ,n k+1 ) (24) 

and 

E{ 1 \Y L ) = E YL (n 1 ,n 2 ) : = 
(Vix? + 1) + V(x? + 2)) + (V(4 1} + 2) + V{x ( l ] + 3)) 

E { k l) {Y L ) = E YL (n k ,n k+1 ) : = 
(V(4 1} + 1) + V{x^ + 2)) + (V(x[ k) + 2) + V(xf + 3)) , (25) 

i.e., these definitions apply to alH = 1, . . . , k. The quantities x^P and (Yl) can be 
thought of as the values of the functions g and E Yl evaluated at the pairs (rii,n i+ i), 
i = 1, . . . , k, of two consecutive gaps in X\, respectively. 

A pair of 2-nd nearest-neighbor molecules in Yl, say the i-th and (i + 2)-th one, is 
separated by two consecutive gaps between 1-st nearest-neighbor molecules, a:^ and a^+V 
From the point of view of the interaction energies between the particles constituting these 
molecules, they are separated by one gap of x\ empty sites, where 

x fi : = xf) + + 2 = {n l/2 } + [n i+2 /2] + (n i+1 + 2). (26) 

Consequently we can define the internal interaction energy of the pair considered of 2- 

(2) 

nd nearest-neighbor molecules in Yl, El (Yl), as the corresponding energy between 1-st 
nearest neighbors, i.e., by replacing simply x± in eq.(P3"D by xf\ 
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While in terms of 1-st nearest-neighbor gaps, rii, the pair of 1-st nearest-neighbor 
molecules, the i-th and (i + l)-th one, was associated with the pair (nj,ni + i), the pair of 
2-nd nearest neighbor molecules, i-th and (i + 2)-th, is associated with two consecutive 
pairs of such pairs ((n*, rij+i), (rij+i, 7^+2)). But from the point of view of the interaction 
energies between the particles constituting the i-th and (« + 2)-th molecules, the transition 
from 1-st nearest-neighbor molecules to 2-nd nearest neighbor molecules can be described 
as a "renormalization process" that maps two consecutive pairs of consecutive 1-st nearest- 
neighbor gaps into a pair of consecutive 2-nd nearest-neighbor gaps: 

((rii, Wi+i), (rii+i, n i+2 )) — ► (rii + (n i+1 + 2), n i+2 + (n i+1 + 2)) . (27) 

Thus the pair of 2-nd nearest-neighbor molecules, the i-th and (i + 2)-th one, can also 
be associated with one pair of gaps - the pair of consecutive 2-nd nearest-neighbor 
gaps {rii + (n i+1 + 2),n i+2 + (n i+ i + 2)). Then in terms of the consecutive 2-nd nearest- 
neighbor gaps 

x f> =x( ni + (rii+i + 2), n l+2 + (n i+1 + 2)) (28) 

and 

Ef\Y L ) = E Yl (m + (ni+i + 2), n l+2 + (n i+1 + 2)) . (29) 

This construction can be naturally continued: the pair of 3-rd nearest-neighbor molecules, 
the i-th and (i + 3)-th one, is associated with the pair of consecutive triples of consecutive 
1-st nearest-neighbor gaps 

((7^, ni+i, n i+2 ) , (n i+ i, n i+2 , n i+3 )) (30) 
or, after renormalizing, with the pair of 3-rd nearest-neighbor gaps 

(m + (ni+i + 2) + (n i+2 + 2), n J+3 + (n i+1 + 2) + (n l+2 + 2)) . (31) 

Consequently 

xf ] =x(m + (n i+1 + 2) + (n t+2 + 2), n l+2 + (n i+1 + 2) + (n l+2 + 2)) (32) 

and 

Ef\Y L ) = E Yl (m + (n i+1 + 2) + (n l+2 + 2), n t+3 + (n l+1 + 2) + (n l+2 + 2)) . (33) 
What we have said above enables us to represent E t ^ lt (Y L ) as follows: 

^r(^) = + E ^ 2) (n) + • • • + s{*>(r L ), (34) 

1=1 1=1 

where the last component in the above equation stands for the internal interaction energy 
of the 1-st and the last, i.e., (k + l)-th molecule in Yj,. 

Having defined a partition of E^iY^) we have to construct a suitable partition of 
E™ l (X), so that the energy differences that will appear eventually in E^iY^) — E™ l (X) 
can be shown to be negative, by means of convexity properties. We have in mind Lemma 
Al of the Appendix. 
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We start with the level of 1-st nearest neighbors. Since in e\ (Yl), given by eq.(|23), 
the energies V(r) group naturally into two pairs, so that each distance involved in one 
pair has a counterpart differing by 1 in the second pair, it would certainly be convenient 
to preserve this property in a partition of E™*(X) we are searching for. Then instead of 
comparing some four energies in E % ^ t (X) with the four energies in E^\Yl) we would have 
to compare only one group of two energies in E^(Yl) with a pair of energies in E l ^ t (X). 

Let us start with the pair V(xf~* + 1), V(x^ + 2) of E^(Yl). In order to be able to 
apply Lemma Al, the corresponding energies in E i ^ t (X) should be associated with some 
distances, say r\ < r 2 , such that r\ + r 2 < (xf^ + 1) + (xf + 2) and r% < xf* + 1 < 
xp + 2 < r 2 . Searching for such distances is facilitated considerably by Lemma A2, which 
gives lower and upper bounds for a quantity like x\. We note that the bounds are the 
best ones, since they can be attained. The particles that constitute the z-th and (i + l)-th 
2-molecules in Yl, whose internal energy is E\ (Yl), come from the z-th, (i + l)-th and 
(z + 2)-th 2-molecules in Xa, separated by the gaps of rij and n i+ i empty sites, or distances 
rij + 1 and + 1, respectively. 

In order not to consider from the very begining the specific behaviour of energies at the 
borders of Xa we limit the index i to the interval 2 < i < k — 1. It will be convenient in the 
sequel to label the particles that constitute the z-th, (i + l)-th and (i + 2)-th 2-molecules 
in Xa by a, b, c, d,e, f, in the order from the left to the right. Thus r a>b = r C)d = r e j = 1, 
where r a >b stands for the distance between the particles labelled a, b, etc. 

In what follows we shall label even (odd) integers by the subscript "even" ("odd"). 
First suppose that rii = n even and n i+ i = n dd- Then by Lemma A4, for n even < n dd 

Ueven + 1 < } + 1 < xf ] + 2 < U odd + 1 (35) 

and 

(n even + 1) + (n odd + 1) = (xf ] + 1) + (x, (1) + 2). (36) 

For n even > n odd , n even has to be interchanged with n odd in ineq. (|35|) . Therefore, according 
to Lemma A4 there is no other choice than n = min{n even + l,n odd + 1} and r 2 = 
max{n even + l,n odd + 1}. Fortunately, there are particles separated by these distances, 
we can set either r x = r& )C and r 2 = r dl e or vice versa. The second group of longer by 1 
distances can be realized as either r x + 1 = r bjd (or r ac ) and r 2 + 1 = r ce (or r d j) or vice 
versa. Now, by Lemma Al 

V(r bc ) + V(r d>e ) = V(n even + 1) + V(n odd + 1) > V(x\ 1] + 1) + V(xf ] + 2) (37) 

and a similar inequality is satisfied if r b c , r d e are replaced by r bid , r c>e (longer by 1) and 
all other distances in ineq. (|37|) are expanded by 1. Summing up, in the case = n even 
and n i+ i = n odd a good candidate for the counterpart of e\ (Y l ) in Xa is 

= (YM + V(r d , e )) + (V(r b4 ) + V(r Cie )) 

= (V(n even + 1) + V(n odd + 1)) + (V(n even + 2) + V(n odd + 2)) . (38) 

Note that only interactions between 1-st nearest-neighbor 2-molecules enter into E^ eo y 
and two one, out of four possible for each pair of 2-molecules, have been picked up. 
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In the second step we suppose that both n» and n i+ i are odd, say n, = n odd and 
n i+ i = m odd . Then by Lemma A4, for n odd < m odd 

nodd + 2 < x^ + 1 < xf ] + 2 < m odd + 1, (39) 

while for n odd < m odd , n odd and m odd have to be interchanged in ineq.(||). If n odd = m odd , 
then 

n odd + 1 = xf ] + 1 < xf ] + 2 = n odd + 2. (40) 
Whatever the relation between n odd and m odd , 

n odd + m odd + 3 = {xf ] + 1) + (xf > + 2). (41) 

Therefore, the distances ri and r 2 are determined uniquely, like in the first step. Trying 
to realize r\ and r 2 as interparticles distances we have to be careful not to run into conflict 
with the choice made previously. The two consecutive gaps of n odd and m odd empty sites 
can be preceded by an even gap, so two interactions between particles belonging to the 
pair of 1-st nearest-neighbor 2-molecules that are separated by the gap of n odd empty 
sites might have already been used. Therefore, in the case Hi = n odd and rii + \ = m odd we 
choose for the counterpart of E^\Yl) in X\ the energy 

Etlo) = (V(n, d ) + V(r dje )) + (V(r a>d ) + V(r c>e )) 

= (V(n odd + 2) + V(m odd + 1)) + (V(n odd + 3) + V(m odd + 2)) . (42) 

In the third step we suppose that 7ij = n odd and n i+ i = n even . Then by Lemma A4, 
for n odd < n even 

nodd + 2 < xr + 1 < x\ 1] + 2 < n even + 2 (43) 

and 

{n odd + 2) + {n even + 2) = {x\ 1} + 1) + (x { - ] + 2). (44) 

For n odd > n even , n even has to be interchanged with n odd in ineq.(|43|). Therefore according 
to Lemma Al there is no other choice than r\ = min{n even + 2, n odd + 2} and r 2 = 
max{n even + 2, n odd + 2}. Therefore in the case Uj = n odd and = n even we choose for 

the counterpart of (Y L ) in the energy 

KL) = (V(n,d) + V(r c , e )) + (V(r a4 ) + V(r cJ )) 

= (V(n odd + 2) + V(n even + 2)) + (V(n odd + 3) + V(n even + 3)) . (45) 

Our definition of Ef) e ^ is not in conflict with those in the previous two situations. 

In the last step we suppose that both and n i+ \ are even, say rii = n even and 
Bj+i = m even . Then by Lemma A4, for n even < m e 



'even 



neven + 2 < xf ] + 1 < xf ] + 2 < m even + 1, (46) 
while for n even > m even , n even and m even have to be interchanged in ineq . (^6|) . If n even = 

neven "I - ! -t~l^Xj ~\~ 2 n even -\- 2. (4^0 



nievem then 
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Whatever the relation between n even and m even , 

Ueven + Eleven + 3 = (xf ] + 1) + (xf ] + 2). (48) 

Therefore, again the distances r\ and r 2 are determined uniquely. As a counterpart of 
E^(Yl) in X\ in the case n« = n e?;en and n i+1 = m e „ en we choose the energy 

^(le) = (V(n, c ) + V{r c>e )) + (K(r M ) + V(r cJ )) 

= (V(n even + 1) + V(m even + 2)) + (V(n e „ en + 2) + V(m even + 3)) . (49) 

This choice is compatible with the fact that the gap preceding n even and the one folowing 
m even can be odd. 

So far we have succeeded in selecting a component of E™ l (X) that is not less than 
X^=2 _1 E\ (Yl)- For any sequence of gaps, n 2 , n 3 , . . . , rik we have assigned, by means of 
definitions (|38"D, (f42|), (f45"|) and fl49|), a sum of pair interactions in such a way that all the 
four pair interactions between the particles constituting 1-st nearest-neighbor 2-molecules 
separated by n i; i = 2,3, ... ,k, empty sites enter the sum. So the sum considered can 
be viewed as the sum over groups of four two-body interactions between particles of 1-st 
nearest neighbor 2-molecules in X. 

It remains to take care of boundary terms. Staying at the level of 1-st nearest neigh- 
bors, in Yl they consist of eight pair interactions: the four ones being the components 
of Ei (Yi) and another four - the components of E^(Yl). The interactions named can 
be compared with the following eight components of E™ l (X): two interactions between 
the left atom and the first 2-molecule, two interactions between the 1-st and the 2-nd 
2-molecules, two interactions between the k-th and (k + l)-th 2-molecules and two in- 
teractions between the k-th 2-molecule and the right atom. It is sufficient to consider in 
detail the left end of A. The interaction of the left atom separated by rti empty sites from 
the 1-st molecule amounts obviously to 

v{n x + 1) + V(m + 2). (50) 

The two still "free" interactions between the 1-st and the 2-nd molecules are: 



V(r bjC ) + V(r a>c ) = V{n 2 + 1) + V{n 2 + 2), if n 2 is odd, 
V(r a:C ) + V(r a4 ) = V(n 2 + 2) + V(n 2 + 3), if n 2 is even. (51) 

Therefore, we define the left boundary internal energy as: 

E { xl = {Vim + 1) + V(n 2 + 1)) + (V(n a + 2) + V{n 2 + 2)) , if n 2 is odd, 

E? e = (Vim + 1) + V(n 2 + 2)) + (V(m + 2) + V(n 2 + 3)) , if n 2 is even. (52) 



Comparing eQ and E\^ with the energies £- i(e/ , r ^, ,,,,„• ^ ,-.„,, , 



,(i) 



E 



(i) 



E 



(i) 



and E, 



H) we find that 



(i) 

i,(e,e) 



defined by 



p(i) 



E 



(i) 

l,(e,o)' 
(1) 

l,(o,o)' 



if ri\ is even 



> E9/„ ^, if n\ is odd 



(53) 
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and 

Thus 
and 



E?l{ = ilf-)' if ^ iseven (54) 

^ > e[% jB) , if mis odd 

E^ - E^fr) > (55) 



EQ-E^(y L )>0. (56) 

At the right end the situation is completly analogous. It is enough to replace n\ by n^, 
^2 by rik + i and interchange the indices e and o in the formulae |53], |54j) to arrive at 
the right boundary internal energies E^l and that are not less than E^(Y L ). 

This lengthy sequence of definitions, given above, amounts to the following: by means 



of definitions (^,0, ^5j^9|,|53i|54D and the remark after inequality (p6|), we defined the 
function Ex on pairs of consecutive gaps, such that its values on pairs of 1-st nearest- 
neighbor gaps, E\ l \x) := Ex(rii,Tii + i), i = 1, . . . , k, satisfy the inequalities 

E\ l \x)-E i i 1 \Y L )>Q. (57) 

Thus J2i=i E^\x) can be taken as the required counterpart of Yli=i E^(Y L ). Now fol- 
lowing the inductive procedure developed to define the partition ( [34]) of E'^ lt (Y L ), we 
apply Ex to pairs of consecutive 2-nd nearest-neighbor gaps and define 

Ef \X) = E x {rii + (n i+ i + 2),n i+2 + (n i+1 + 2)) (58) 

By Lemma Al the inequality (f>7\ ) implies that 

EP(X)-EP(Y L )>0, (59) 

so Yli=i (X) can stand for the required counterpart of VjlZ\ E^ (Y L ) . Continuing this 
process we arrive at the following partition of E l ^ t {X): 

k k—1 / fc+1 \ 

Ef\X) = £ E<?\X) + Y: E?\X) + ... + E{ k) (X) + W m + + 2) + 1 (60) 

i=l i=l \ i=2 I 

and moreover we obtain that 

ET\X) - Et\Y L ) >v(nx + + 2) + 1^ > 0, (61) 

where V + 2f=2 ( n i + 2) + lj stands for the interaction of the left atom with the right 
one. Clearly, if we substitute Yr for Yi we arrive at the same conclusion as fl6~I|) . 

Remark: The inequality ( Ell ) holds true if there are no molecules between two atoms. 
Lemma 2 (Eliminating of n-molecules with n > 3) 

Consider a lattice gas ([[]) with an interaction energy V 2 , such that V 2 (l) > 0. Let X 
be a configuration containing triples of particles separated by distance one. Then there 
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is a configuration Y, with p(Y) = p(X), free of such triples. Moreover, if V 2 (2) > 
(7W - V 2 (l))/2, then e(Y) < e(X), where I¥ = E^ 3 V 2 {r). 

Corollary: In a lattice gas of Lemma 2, if V 2 (2) > (7W — V 2 (l))/2, then the ground-state 
configurations do not contain triples of particles separated by distance one. In particular, 
the ground-state configurations do not contain n- molecules with n > 3. 

Proof: The idea of our proof is to start with an arbitrary configuration containing forbid- 
den triples of particles, that is triples of particles separated by distance one, and transform 
it, using Lemma 1 and some rearrangements of particles, whose energy costs are controled 
by simple bounds, into another configuration of the same density but free of the forbid- 
den triples of particles. The proof is completed by showing that the overall energy cost 
of attaining the final configuration is negative, provided that a suitable condition on the 
interaction energy is imposed. 

Let X be an arbitrary configuration with the particle density p(X) < 1/2 and whose 
density of the forbidden triples, P3(X), satisfies the inequality pz{X) > 0. Removing one 
particle from each forbidden triple we lower the energy by at least V 2 (2) + V 2 (l) and arrive 
at a configuration X whose particle density satisfies the double inequality 

l/2>p(X)>p(X)-p 3 (X). (62) 

Consider the class of configurations that consist of 2-molecules and atoms (separated by 
at least two empty sites) exclusively and whose density is p(X). The configuration X 
belongs to this class. Let Y be the lowest energy density configuration in this class, thus 
e(Y) < e(X). By (p2p, p(X) — p(X) < p 3 (X), therefore the energy density variation 
e(X) - e(X) is bounded from below by (V 2 (2) + V 2 (l)) (p(X) - p(X)). Consequently, we 
arrive at the following lower bound for e(X): 

e(X) > (V 2 (2) + V 2 (l)) (p(X) - p{X)) + e(Y). (63) 

By Lemma 1 the configuration Y does not contain atoms, thus it consists exclusively 
of 2-molecules. It is clear that there is only one (up to translations) configuration of 
density 1/2 that consists exclusively of 2-molecules. It is the periodic configuration whose 
elementary cell has the form [• • oo], where the filled circles stand, say, for occupied sites. 
The 2-molecules are separated by distance 3, i.e., by 2 unoccupied sites. Obviously, a con- 
figuration that consists exclusively of 2-molecules and whose particle density is less than 
1/2, contains arrangements of 3 unoccupied sites in a row with a nonvanishing density. 
Such a configuration can be thought of as a collection of local particle configurations that 
are separated by at least 3 unoccupied sites in a row. Shifting the local configurations 
separated by at least 3 unoccupied sites in a row, while respecting the conditions that 
they cannot be separated by less than 2 unoccupied sites in a row, we can create arrange- 
ments of 6 unoccupied sites in a row. Then in the middle of the 6-site gap obtained we 
insert a 2-molecule. The whole procedure is repeated untill the particle density comes 
back to the value p(X). The result is a configuration Y. Now we can try to compare the 
energy densities e(Y) and e(X) and see if e(X) can be strictly larger than e(Y), what 
would prove that ground-state configurations do not contain forbidden arrangements of 
particles. Suppose that we can bound from above by A the energy cost of creating a 6-site 
gap and inserting there a 2-molecule. Then we can compare e(Y) and e(Y): 

e(Y)-e(Y)<A(p(X)-p(X))/2. (64) 



19 



By means of inequalities (|62| ) and (|63"D, we arrive at the folowing inequality relating e(X) 
and e(Y): 

e(X) > e(Y) + (V 2 (2) + V 2 {1) - A/2)(p(X) - p(X)). (65) 



Therefore, the ground-state configurations do not contain forbidden arrangements of par- 
ticles if V 2 (2) + Va(l) > A/2. 

It remains to derive the upper bound A. For that we first estimate the energy cost of 
translating a single particle by distance 1. Let a chosen particle be separated from its left 
nearest neighbor by distance r and from its right nearest neighbor by distance r . The 
total interaction energy of the particle can be viewed as the sum of the interaction energies 
E^\r) and E^'(r) with all the particles to its left and to its right, respectively. Suppose 
we shift our particle by distance 1 to the left. Then the change in the total energy of the 
particle is bounded from above by E^\r — 1) — E^\r), since E^{r + 1) — E^\r') < 0. 
Let Xi, i — 1, 2, . . ., be the distance between the i-th left nearest neighbor and the (i + l)-th 
left nearest neighbor. Then 

E%\r - 1) - E { l\r) = (V 2 (r - 1) + V 2 (r - 1 + x x ) + V 2 (r - 1 + x x + x 2 ) + . . .) 

- (V 2 (r) + V 2 (r + xi) + V 2 {r + x x + x 2 ) + . . .) = V 2 (r - 1) + 
- (V 2 (r) - V 2 (r - 1 + xi)) - (V 2 (r + x x ) - V 2 (r - 1 + x x + x 2 )) -...<V 2 (r- 1). (66) 

By means of inequality ( p6| ) we are able to give a simple upper bound for the energy 
cost of translating a local configuration considered above by distance 1 to the left (right). 
This energy does not exceed the sum of V 2 (r — 1) over all distances r, r > 4, between 
the particles of the configuration considered and the left nearest neigbor particle of this 
configuration: 

oo oo 

n r V 2 (r - 1) < E V 2 (r - 1) = W, (67) 

where n r = 1 if the local configuration considered contains a particle separated by distance 
r from the left nearest-neighbor particle of this configuration, and n r = otherwise. 
Consider again a configuration that is a collection of local configurations separated by 
at least 3 unoccupied sites. We can create a 6-site gap between two local configurations 
by shifting a chosen local configuration by distance one to the left, then shifting its right 
nearest-neighbor local configuration by distance 1 to the right, what creates a 5-site gap 
between the configuration chosen and its right nearest-neighbor local configuration, and 
after that shifting simultaneously the right nearest-neighbor and the right next nearest- 
neighbor local configurations by distance 1 to the right. Each of the three reshuffles of 
local configurations increases the energy by no more than W. Thus, a 6-site gap can be 
created at the cost that does not exceed 3W. Now inserting a 2-molecule in the middle 
of the 6-site gap created will cost no more than AW + ^(1). Therefore, we can set 
A = 7W + V 2 (1). 
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5 Some extensions of the main result and applica- 
tions 



The both lemmas, formulated and proved in the previous section, set conditions on the 
value of Vi at distance one, which go in opposite directions. It is for this reason that in 
Theorem 1 this value is set to zero. However, it is not hard to see that with some effort 
this condition can be relaxed. 

First, suppose V 2 (l) < 0. This condition fits the assumptions of Lemma 1 but not 
those of Lemma 2. In the latter case, the energy gain on breaking forbidden triples of 
particles is as before V^(2) + V 2 (l) but is diminished now by the second term and may even 
become negative. On the other hand, the energy increase on inserting a 2-molecule into a 
six-site gap, 4VF+V2(1), is also diminished. Thus, Theorem 1 holds true also for V 2 (l) < 0, 
but such that V 2 {2) > (7W - V 2 (l))/2, which can only take place if |V^(1)| < 2V 2 (2). 

Second, suppose V 2 (l) > 0. Now, in turn, this condition fits the assumptions of 
Lemma 2 but not those of Lemma 1, which is a more severe problem to overcome than 
the previous one. The final conclusion of our proof of Lemma 1 is a consequence of the 
fact that the energy difference of an initial local configuration, consisting of two atoms 
separated by some k 2-molecules, and the final one, consisting of k + 1 2-molecules, is 
not smaller than V 2 (r atom ) — V 2 (l), where r atom is an arbitrary distance between the two 
atoms in the initial local configuration. If, as we have supposed originally, V 2 (l) = 0, then 
the lower bound obtained implies that the energy difference considered is positive what 
completes our proof of Lemma 1. If however V 2 (l) > 0, to arrive at the same conclusion 
by means of a similar proof we would have to limit from above r atom , which is not possible. 
Instead, we can turn to a strategy analogous to that used in Lemma 2. 

In order to restrict the set of possible ground-state configurations to the set C\ p we 
have to get rid of forbidden pairs of particles, i.e., particles separated by distance 1. 
We fix distance r > 4 and choose a particle density p < 1/r. Removing one particle 
from each forbidden pair we lower the energy by Vi(l) and arrive at a configuration 
of lower density, consisting of atoms exclusively. The ground-state configuration among 
such configurations is a most homogeneous configuration, in which the smallest distance d 
between neighboring atoms is greater than r. Then we insert back the removed particles in 
such a way that the smallest distance between neighboring atoms is greater than [d/2] + 1. 
Now, if Vi(l) > 2 J2 r >[d/2]+i ^i( r )> then the obtained configuration has a lower energy 
density. Therefore, under the stated conditions, any ground-state configuration consists 
only of atoms and consequently is a most homogeneous configuration. 

Finally, we would like to point out that our Theorem 1 and its generalizations discussed 
above can be used to study properties of two- and three-dimensional lattice gases. 

Consider first a two-dimensional lattice gas on the square lattice Z 2 . Let (x,y) be 
an orthogonal coordinate system, such that the its axes coincide with lattice directions 
with unit lattice constants. Suppose that with respect to this coordinate system the 
Hamiltonian decomposes into H\ + H\, where is a sum of two-body potentials 
supported by pairs of sites along the re-direction, and H\ is an analogous sum but in 
the ^/-direction. Then, knowing the one-dimensional ground-state configurations of H^, 
we can easily construct two-dimensional ground-state configurations of H\. For instance, 
the particle configurations of lattice lines in the x-direction can be chosen as the ground- 
state configurations of H^, with the condition that, going in the positive y-direction, the 
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configuration of the next lattice line is shifted by one lattice constant in the positive 
x-direction. The Hamiltonians considered above naturally appear in surface physics [TIJ . 



Another situation, where our results are applicable, is the realm of three-dimensional 
layered systems, like those studied extensively by Fisher and collaborators, see for instance 
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6 Appendix 

Definition 1 We say that a real function /, defined on integers n > no, is convex, if for 
every n > n + 1 the following inequality is satisfied: 

f(n + l)+f(n-l)>2f(n). (68) 

Lemma Al (Convexity of decreasing / in Z) 
Let integers x, y, s, t satisfy the inequalities 

Uq< x < s <t <y and x + y < s + t (69) 

and let / be a convex and decreasing function for integers n > n a . Then 

/(*) +/(») >/(*)+ /(*)■ (70) 

Definition 2 For any positive integer n, [n/2] is the greatest integer that does not exceed 
n/2 while {n/2} = n - [n/2]. 

Lemma A2 (Properties of bracket functions [.] and {.}) 
For positive integers n < m, [|] < [y], and {|} < {y}. 
Moreover, 



m + n 



2 



[y] + [|] + 1, if m and n are odd 



. 2 J 1 L 2 J 

"' 1 " 1 otherwise 



and 

m + n\ _( {y} + {2} — 1 ; if m and n are odd 
I {f } + il }> otherwise. 



(71) 
(72) 



Definition 3 For an ordered pair of positive integers (s, t) we define 

g(s,t) = {s/2} + [t/2]. (73) 

Lemma A3 (Properties of function g) 
For positive integers s <t, 

s<g( s ,t)<t. (74) 
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More refined bounds are given in Lemma A4. For positive integers p, s, t 

g(s+p,t+p)=g(s,t)+p. (75) 
Lemma A4 (Best bounds for function g) 

Let n even and m even be two positive even integers and n odd and m odd - two positive and 
odd integers and let g e>0 be an abreviation for g(n even ,n dd), etc. Then 

EO : (n even , n od d) 

(EOl) If n even < n dd, then 

Tleven — 9e,o ^ f^odd 1 

(E02) If n even > n odd , then n odd < g e ,o < n even - 1 

2g e ,o = n even + n odd - 1 (76) 

00 : (n odd , m odd ) 

(001) If n odd < m odd , then n odd + 1 < g a ,o < m odd - 1 

(002) If n odd > rn odd , then m odd + 1 < g ,o < n odd - 1 

(003) If n odd = m odd , then g 0>0 = n odd = m odd 

2g ,o = m odd + n odd (77) 

OE : (n odd , n even ) 

(OEl) If n odd < n even , then n odd + 1 — 9o,e — ^even 

(0E2) If n odd > n even , then n even + 1 < g 0>e < n odd 

2g , e = n odd + n even + 1 (78) 

EE . {jleveni Eleven) 

\EE\j If n even <c m eveny then n even -\- 1 ^ 9e,e — vrieven 1 

(EEQi) If n even > m even) then ni even -\- 1 ^ 9e,e — ^even 1 

2f/ , e,e THeven ~\~ IT'even (^) 
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